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We introduce a method of intervals for the analysis of diffusion-limited annihilation, A + A — > 0, 
on the line. The method leads to manageable diffusion equations whose interpretation is intuitively 
clear. As an example, we treat the following cases: (a) annihilation in the infinite line and in infinite 
(discrete) chains; (b) annihilation with input of single particles, adjacent particle pairs, and particle 
pairs separated by a given distance; (c) annihilation, A + A — > 0, along with the birth reaction 
A — > 3^4, on finite rings, with and without diffusion. 
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I. INTRODUCTION 



Diffusion- limited reactions display a wide range of behavior characteristic of noncquilibrium dynamics, such as 
self-organized criticality, pattern formation, and dynamic phase transitions. However, the exact description of this 
behavior can only be determined for the simplest reaction schemes. Single-species annihilation, A + A — ► 0, and 
coagulation, A + A — > A, are among the few examples of exactly solvable diffusion-limited reactions. The Method 
of Empty Intervals has provided many exact results concerning diffusion- limited coagulation . Its advantage lies 
in the fact that it leads to simple diffusion equations which are easy to solve, and whose interpretation is intuitively 
obvious. We develop a similar method of intervals for the study of diffusion-limited annihilation reactions. 

Diffusion- limited annihilation has been studied in a variety of forms. For the basic annihilation process, the exact 
time-dependent concentration of particles has been determined for various initial distributions on the infinite one- 
dimensional lattice and the continuous real line p|-[il| . The inclusion of particle input in the reaction scheme allows for 
a nontrivial steady state concentration. Both input of adjacent particle pairs and single particle input |p7| , ^9| , p0| 

have been studied. A nontrivial steady state may also occur when a birth reaction is included. Annihilation with 
the symmetric birth reaction A — > 3A has been studied for its implications to interacting particle systems and 
universality theories [^l|-^8|. Sudbury ]23j referred to this reaction as the Double Branching Annihilating Random 
Walk (DBARW). Without diffusion, the process is known as the Double Branching Annihilating Process (DBAP). 
We show how a method of intervals, previously introduced in the literature |^ , p9| -p2| , may be used to model all of 
the above reactions. This powerful yet simple method yields insights that extend beyond the known results p3[ . 

The rest of this article is organized as follows. In section [n] we introduce the method of intervals, as adapted for 
the annihilation process A + A — > 0, and reproduce its well-known kinetics. Section [II deals with various cases of 
input. Input of particle pairs separated by a distance y is considered, and the limiting cases of input of single particles 
(y —>■ oo) and adjacent particle pairs (y — > 0) are examined. We also discuss input at a rate proportional to a function 
of the global concentration of particles. In section IV we present our results for DBARW and DBAP. We conclude 
with a discussion and open directions of research suggested by the new method of intervals, in section [v| 



II. THE METHOD OF ODD/EVEN INTERVALS 



We now adapt the method of empty intervals, conceived originally for the analysis of diffusion-limited coalescence, 
A + A — > A, to reactions such as annihilation, A + A — > 0, where the parity (the number of particles modulo 2) is 
conserved. Consider the process A + A — > 0, taking place in a one-dimensional lattice. The particles hop randomly 
to the nearest site on their right or left, at equal rate T, and annihilate immediately upon encounter. Let G n (t) be 
the probability that an arbitrary segment of n consecutive sites contains an even number of particles, at time t. (We 
assume, for the moment, that the system is infinite and homogeneous.) A site can be either empty or occupied by a 
single particle, so the probability that a site is occupied, i.e., the particle density, is 

P (t) = l-G 1 (t) . (1) 
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Since the reaction A + A — > conserves parity, the only way that G„ might change is when particles at the edge 
of the segment hop outside, or when particles just outside of the segment hop inside. To describe these events, we 
require F n (t) — the probability that an n-segment containing an even number of particles is followed by the presence 
of a particle at the (n + l)-th site. This can be expressed in terms of the G„ (Figure 1): 

2F n (t) = (1 - Gi) + (G n - G n+ i) , (2) 

Likewise, H n (t) — the probability that an n-segment containing an odd number of particles is followed by a particle 
at the (n + l)-th site — is 

2H n (t) = (1 - G x ) - {G n - G n+ i) . (3) 
The evolution equation for G n is now readily obtained: 

^-G n (t) = 2T(F„_ 1 - + H n - F n ) . (4) 

at 

The term proportional to F n ^\ on the r.h.s. expresses the event that a particle at site n jumps out of the segment, 
leaving an even number of particles in the remaining n — 1 sites (and hence in the n-segment); H n -\ corresponds to 
the same case, but when there arc initially an even number of particles in the n-segment (that is, an odd number in 
the (n — l)-segment); H n and F n pertain to a particle just outside of the n-segment, at site n + 1, jumping in. The 
factor of 2 accounts for events taking place at both edges of the segment, at equal rate. Using Eqs. (^) and (|J), this 
becomes 

jG n {t) = 2r(G„_ 1 -2G„ + G„ +1 ) . (5) 

The case of n = 1 requires a special equation, since Go is undefined. Taking into account all the ways Gi might 
change, one finds 

^Gi(i) = 2r(l-2G! + G 2 ) . (6) 
ot 

Thus, Eq. (^) may be understood to be valid for all n > 1, provided that one uses the boundary condition 

G (<) = 1 • (7) 

Additionally, since the G„ are probabilities, we have 

< G n (t) < 1 . (8) 

The initial distribution of particles determines the required initial condition, G n (0). For example, suppose that at 
the start of the process each site is occupied randomly and independently of other sites, with probability pq. Then 

Gn(0) = i + i(l-2 Po r. (9) 

Eq. (|^) may be solved by standard techniques, for example, by Laplace-transforming with respect to time, fitting an 
exponential solution to the resulting difference equation, and finally inverting the Laplace-transformed solution p4[ . 
With the boundary conditions (Q) and (||), and the natural initial condition (||), one obtains: 

G n (r) = \ + \(1 ~ 2po)"e^ + | jf e^' (1 - e^')I n (2r') *±- , (10) 

where r = 2rt, (3 — 4/5q/(1 — 2po), and /„(•) is the modified Bessel function of order n |3j|. In particular, the 
probability that a site contains a particle, is 

p(r) = 1 - G x {r) = \ - \(1 - 2p )e&- -l£ e^' (1 - e^')h(2r') ^ . (11) 

It is often more convenient to assume a low initial density of particles, and work in the continuum limit. Let a be 
the lattice spacing. Then, putting x = na, G n (t) — > G(x,t), and T = D/a 2 in Eq. (|^), and taking the limit a — * 0, 
one obtains the diffusion equation 

-G(x,t) = 2D—G(x,t) , (12) 



2 



with the boundary conditions 



The particle concentration is obtained from 



G(x,0) = l, (13a) 
0<G(x,t)<l. (13b) 



c(t) = lim^ = -^-G(a; ) t)| s= o. (14) 

Consider, for example, an initial concentration cq of randomly placed particles. The initial condition is obtained 
from (|^), putting p Q = c Q a and passing to the continuum limit: 

G{x,Q) = \ + \e-^ . (15) 
Solving for the concentration is then straightforward: 

c(t) = c e z2 erfc(z) , z = 2c V2Dt . (16) 
This is the known result, with its familiar long-time asymptotic behavior, 

ctt) ~7mi- (17) 

III. ANNIHILATION WITH INPUT 

Let us now include the possibility of input, where empty sites become spontaneously occupied at a prescribed rate. 
We shall assume that the input of particles is homogeneous (translation invariant). Various cases of input may be 
analyzed through the method of intervals within this restriction. 

A. Input of single particles 

Consider a random, homogeneous input of single particles at constant rate R per unit space per unit time. If the 
input site is already occupied we assume that the particles react immediately and the site becomes empty. In other 
words, the state of individual sites flips (from empty to occupied, and vise versa) at rate Ra per unit time. This kind 
of input affects the rate of change of G n thus: 

(^-) = nRa(l - G n ) - nRaG n = nRa(l - 2G n ) . (18) 

V Ot I A-input 

The term nRa(l — G n ) accounts for the increase in G n due to the input of a particle to an n-site interval, initially 
containing an odd number of particles, while ~nRaG n is the (negative) change in G n when a particle is input into 
an initially even n-interval. Note that Eq 
condition (||). In the continuum limit, the c 

—G(x, t) = 2D—G(x, t) + Rx(l - 2G(x, t)) , (19) 

subject to the same boundary conditions as without input. The problem can be solved by substituting 

K(x,t) = 1 - 2G(x,t) . 

K(x, t) then satisfies 

j t K{x, t) = 2D^K{x, t) - 2RxK(x, t) , 
K(0,t) = -1 . 



lq ) is valid for all n > 1, and so the input does not affect the boundary 
lange due to input is dG/dt = Rx(l — 2G), so Eq. (Oh becomes 
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Using the method of separation of variables, we write 

A>0 



where the eigenfunctions satisfy 



Thus, 



2D-^K x {x) = {2Rx-\)K x (x) 



Kx(x)=M 



i?\l/3 A 



D) ' 2(DR 2 y/ 3 



where Ai(-) is the Airy function — a solution to the equation Ai"(z) = zAi(z) |3f|. The boundary condition K(0, t) = 
— 1 implies Kq(0) = — 1 and K\(0) = (A > 0). The condition A'o(0) = — 1 yields the steady state solution: 

G ^-2 + 2Ai(0) ' (20) 



d rM , 1^(0)1 ^1/3 
c s = —£-G a {x)\ x=0 = - - — . (21) 



and the steady state concentration: 

The condition ^a(O) = imposes discrete eigenvalues: 

A„ = 2( J Di? 2 ) 1 / 3 |a„| , 

where a n are the zeros of the Airy function. The relaxation time to the stationary state, r, is then given by the 
minimal eigenvalue: 

r" 1 = A min = 2(L>i? 2 ) 1 / 3 |ai| = 4.6762(Di? 2 ) 1 / 3 . (22) 

B. Input of adjacent particle pairs 

Consider now homogeneous, random input of particle pairs to adjacent sites: any two adjacent sites become occupied 
at rate Ra = r per unit time. As before, if a target site is already occupied, the site becomes empty as a result of 
input. The situation is analogous to that of input of single particles, but the parity of an interval is now only affected 
when input occurs at the interval's edge — when only one particle of the pair falls right inside the interval. G n is 
increased by input at the edge of an odd n-interval, and decreased by input at an even interval (See Figure 2). Thus, 
the rate of change of G n due to pair input is 

(^-G n ) = 2[r(l - G n ) - rG n ] = 2r(l - 2G n ) . (23) 

V Ot I AA-iuput 

In the continuum limit, G(x,t) must satisfy 

-G(x, t) = 2D—G(x, t) + 2r(l - 2G(x, t)) , (24) 

with the usual boundary condition, G(0, t) = 1. 

Again, we use K(x,t) = 1 — 2G(x,t) to determine the solution. K(x,t) must satisfy 

^K(x, t) = 2D^K(x, t) - ArK(x, t) , 
K(Q,t) = -1 . 
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Expanding K(x,t) — Y^\K\(x)e , as before, we find 

2D-^K x {x) = (4,r-X)K x {x) , 
with K (Q) = — 1 and K\(Q) = (A > 0). This leads to the stationary solution 



and the stationary concentration 



2 \D 

The relaxation spectrum is continuous: 



G.(X) = j + , (25) 



A > 4r ; = sin(^/ ^^^) , (27) 

and the relaxation time is 

r- 1 = A min = 4r . (28) 

These results, along with the single particle input results, are in complete agreement with previous work by Racz [T^ ] 
and Lushnikov [ITsfl . 



C. Input of correlated pairs 



The method of intervals allows us to deal with more complicated input. As an example, we consider input of particle 
pairs, separated by to lattice spacings. This input of correlated pairs i nterpolates between the t wo cas es discussed so 
far: for m = the particle pairs are adjacent, just as in section IIIB. while the case of section III A is recovered as 
to — > oo, since then the input particles cannot possibly affect each other and the correlation is lost. 

Suppose that the particle pairs are deposited at rate Ra per unit time (per site). The rate of change of G n due to 
input is now different for n < m and n > to: 



(!°«L 



A-input 



A-input 



2Rna(l - 2G„) , n < to , 

2R(m + l)a(l - 2G n ) , n>m. 



Again, Eq. (|29|a) does not affect the usual boundary condition Gq = 1. Adding these contributions to Eq. 
taking the continuum limit, we obtain 



^G(x, t) = 2D^G{ Xl t) + 2Rx(l - 2G(x, i)) , 



dt 




dx 2 
d 2 



gi G(x, t) = 2D—G(x, t) + 2Ry{\ - 2G(x, t)) 



x < y , 



x > y , 



with the usual boundary condition G(0,i) = 1. Additionally, G{x,i) and dG{x,t) / dx are continuous at x = y. 
For simplicity, we analyze only the steady state. We find, 



G s (x) = c%Ai(ax) + C2Bi(ax) 
G s (x) - c 3 e-V^* + I , 



1 



x < y , 

x > y 



(29a) 
(29b) 
]) and 

(30a) 
(30b) 

(31a) 
(31b) 



where Bi(-) is the independent, divergent solution to Airy's equation |3j|, a = (2i?/I?) 1 / 3 , and ci, C2, C3, are constants 
obtained from the boundary conditions at x — and x — y. The steady state concentration is 



_ lAi'(0)| /2flxV3 
s ~ 2Ai(0) V D ) 1 



y 



(32) 
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where 



^Bi(z) + Bi'(z) + VTzAi(z) + V3Ai'(z) f^, z « 1 , 
V5Bi(z) + Bi'(z)- VSzM(z)- V3Ai'(z) ^ \l , z»l. 

Eq. ( |33| ) expresses the crossover behavior between input of single particles and adjacent particle pairs. Indeed, when 
y — > oo (z — > oo), the particles are essentially uncorrelated, and the steady state is the same as for single particles, of 
Eq. , but with 2R instead of R. This is because the input of pairs introduces particles at twice the rate of single 
particles. The limit of y — > (z — * 0) yields c s ~ \J (R/D)y, just as for adjacent particle pairs (Eq. ^). The crossover 
between the two regimes occurs about (R/D^^y ~ 1 (Figure 3). 



D. Input proportional to some global property 

We can also analyze inpu t at a rate proportional to some global property of the system. Consider input of single 
particles, as in section [II A, but at a rate proportional to a functional of G(x,t): 

R = n[G(x,t)] . 

Putting this rate into Eq. ( p~9| ) yields a non-local partial differential equation for G(x,t), which is generally difficult 
to solve. However, the steady state may be obtained in the following manner. In the steady state limit, the rate 
R s = TZ[G s (x)] is constant. This constant can be found by solving 

R s = TZ[G s (x\R s )}, (34) 

where G s (x\R) is the distribution given by Eq. (p0|), assuming constant rate R. Once the value of R s is known, the 
steady state distribution, G s (x\R s ), and steady state concentration, c s , follow. Notice that the same procedure applies 
for other kinds of input, with minor changes: for example, for input of adjacent particle pairs, G s (x\R) is given by 
Eq. (H) rather than @. 

An amusing example is input of single particles at a rate proportional to a power of the particle concentration: 

* = <»> 
In this case we can find c s directly, from Eq. (pl|): 

_ |Ai'(0)| 

7 2Ai(0) ' 



\-d) ' 



or, 



l^^J • (36) 



,c(0) 

The true meaning of Eq. (|36|), for different values of a, is revealed only from a careful analysis of the pertinent 
effective-rate equation Q|: 

j f c = -k lC 3 + k 2 c a , fei ~ D , k 2 ~ R o /c{0) a . (37) 

Eq. (p7|) has a steady state similar to (^), but one can clearly see that that state is stable only for a < 3. For 
a > 3 the solution is unstable — instead, the system flows to one of the stable states c = or c = oo. A particularly 
intriguing case, in which fluctuations not reflected in the mean-field analysis will surely play a central role, is that of 
a = 3. 



IV. ANNIHILATION WITH SYMMETRIC BIRTH 



We now treat the process of annihilation, A + A — > 0, with the addition of the back reaction A — > 3A, where particles 
give birth to two new particles at the adjacent sites to their left and right, at rate f2 (per particle, per unit time). If 



G 



birth takes place onto a site that is already occupied, annihilation is immediate and the site becomes empty. This 
process has been analyzed by Sudbury |^| , who named it the double branching annihilating random walk (DB ARW) . 
Sudbury pj| also studied the double branching annihilating process (DBAP), which is the same as the DBARW except 
that the particles do not diffuse. The DBARW is now commonly referred to as the branching-annihilating walk with 
two offspring (n — 2-BAW), but since we wish to distinguish between the processes with and without diffusion we 
shall utilize Sudbury's nomenclature. 

Consider the DBARW taking place on an iV-site chain, with periodic boundary conditions (or an iV-site ring). 
Because the ring is finite, Gat is determined by the parity of the initial number of particles, and remains constant 
throughout the process. The effect of diffusion on G n has already been discussed, and is given by Eqs. (§) and (|): 

(^G n ) = 2T(G n . 1 -2G n + G n+1 ) , Kn<N, 

\Ot J Diff 



and 

d 



(-gA =2r(l-2Gi + G 2 ) , n = l 

\0t ) Diff 



Birth affects G n in a similar way to diffusion: the parity of an n-segment changes only when a particle just inside or 
just outside the segment gives birth. Thus, 

( ^-G n ) = 2ft(G„_! - 2G„ + G n+1 ) , K n < N - 1 . 

\Ot / Birth 

For n = 1, birth into an empty site decreases G\, while birth into an occupied site increases G\. The two effects add 
up to 

d 



(ttGi) = 2ft(G 2 - Gi) , n = 1 

V ut J Birth 



Finally, for n = N — 1, birth from a particle at the inner edge of the (N — l)-segment is similar to the case of generic 
n, but birth from a particle outside of the segment is different: because the lattice is a iV-ring, a particle outside of 
the (iV — l)-segment gives birth to two particles inside it (one at each edge), and the parity does not change. Thus, 



( — Gjv-i) — 2CI(Gn-2 — Gn-i) 

V ot J Birth 



n = N - 1 



Putting the contributions from diffusion and from birth together, we get 



^G„=2(r + 0)(G„_ 1 -2G„ + G n+1 ), Kn<A-l, (38a) 

-^Gi=2r(l-2Gi + G 2 ) + 2fl(G 2 -G 1 ) , n = l, (38b) 

^G A r_ 1 =2r(G A r_ 2 -2G A r_ 1 +G A r)+2r!(G A r_ 2 -GAr_ 1 ) , u = N - 1 , (38c) 

For simplicity, we focus on the steady state, where time derivatives are zero. (The transient behavior may be 
analyzed by methods similar to that of Section ||.) The general (steady state) solution of Eq. (|38|a) is G„ = An + B, 
where A and B are constants. Their explicit value is determined from the boundary conditions (|38|b) and (|38|c): 

Q.A - TB = -r , (39a) 
[NT + Q)A + TB = -TG N . (39b) 

Thus, 

A= F(1-GaQ ft(l-GW) 
AT + 2n ' ' AT + 2fi ' 

and the stationary particle density is 



7 



When the initial number of particles is even, Gn = 1 and the system gravitates into its absorbing, empty state — a 
state from which it cannot evolve any further. If the initial number of particles is odd (Gn = 0), the system can never 
reach the absorbing state. However, in this case the steady state density is barely larger than 1/N (for N large): the 
system comes as close to extinction as possible. This generalizes the result of Sudbury that p s — for the DBARW 
in infinite lattices. 

Consider now the DBAP, which is the case of no diffusion, T = 0. In this case, Eqs. (|39| ) imply A = 0, B arbitrary, 
or G n = const. Suppose that the initial number of particles is infinite, so that Gn is not defined (but, since G n is 
constant, < Gn < 1)- Then, taking the limit T — > in (ff0|), we find < p s < 1/2. This too agrees with Sudbury, 
who has shown that any homogeneous, random initial distribution leads to p s = 1/2, but other distributions lead to 
steady states with < p s < 1/2. According to Sudbury, the possible steady states are the borders of the product 
measures v p . Indeed, the border of v v has constant G n : the condition that there be an even number of borders in an 
n-segment is equivalent to having the same states at the edges of the segment in the original v v measure. But the 
probability of this event is G n = p 2 + (1 — p) 2 , independent of n. 

An amusing case of the DBAP is when the initial number of particles is finite. In this case Gn is known exactly. 
Thus, if the initial number of particles is odd (Gn — 0), then the r — + limit of Eq. ( ^p|) yields p s = 1/2, in agreement 
with Sudbury. It is instructive to see how this steady state manifests itself in the case of finite rings. An example, of 
N = 6, is shown in Appendix A. If the initial number of particles is even (Gn = 1), Eq. ( [lo|) suggests that p s = 0, 
even for finite N. This is however not true, and it can be shown p^ , |2^ ] that the number of particles remains bounded. 
(For example, two adjacent particles would propagate as a pair forever, diffusing without change.) In this case the 
limit r — > is singular: even the tiniest amount of diffusion would land the system in the absorbing state, but the 
system can never become empty without diffusion. 



V. DISCUSSION 



In conclusion, the method of intervals enables one to obtain exact results for a large class of diffusion-limited 
annihilation models in one-dimension. The usefulness of this approach stems from the direct consideration of parity 
conservation in the annihilation process; annihilation conserves the number of particles modulo 2. Defining G n as 
the probability that n consecutive sites contain an even number of particles exploits this constraint. In some sense, 
this method of intervals is a generalization of the method of empty intervals used for coagulation processes [fl| — |5f1 - A 
further generalization exists that simultaneously models the g-state Potts model, coagulation and annihilation [ |3l[ . 

The results of this paper can be summarized as follows. The well known result c(t) = Co e 8c ° Dt erfc ( y8c| Dt) for 
diffusion-limited annihilation with an initial Poisson distribution of particles in the real line was reproduced, and we 
have obtained an analogous solution for the discrete case of a linear lattice. We have studied the nontrivial steady 
state concentration resulting from particle input. In particular, it was shown that a smooth crossover occurs in the 
steady state concentration as the separation of input particle pairs increases from zero to infinity; as particle pair 
input crosses over to single particle input, the concentration dependence changes from (R/D) 1 ^ 3 to (R/D) 1 / 2 . Our 
final example involved annihilation with the symmetric birth reaction A — ► 3A The method of intervals allows one 
to fully examine the kinetics of this reaction scheme in finite ringed lattices. As an example, we derived the exact 
nature of the steady state solution for even and odd populations, with or without diffusion. 

The method of odd/even intervals introduced here for the case of annihilation could be extended along the same 
lines as the method of empty intervals used for coagulation j^||||56],[37|] . This will be the subject of future work. 
For example, the kinetic phase transition associated with reversible coagulation, A + A ^ A pq] , might also have 
an echo in A + A — > 0, A — + 3A, at least in the case of finite lattices, when there exists a non-empty steady state. 
Nonhomogeneous systems could be studied with the addition of one variable, by focusing on G nm (t) — the probability 
that the interval between sites n and m contains an even number of particles . It might also be possible to obtain 
multiple-point correlation functions by studying the joint probability that several distinct intervals contain odd/even 
numbers of particles J|||. 

A natural question to ask is whether the inter-particle distribution function (IPDF) for annihilation can be computed 
by the odd/even interval method, as is the case for coagulation (with the method of empty intervals). In the 
latter case, the IPDF — the probability that the empty space between two particles is of length x — is given by 
p(x, t) = c(t)~ 1 d 2 E(x, t)/dx 2 , where E(x, t) is the probability that an interval of length x is empty Q. Unfortunately, 
d 2 G(x, t) I dx 2 does not convey analogous information. We believe that the full hierarchy of multiple-point correlation 
functions may be obtained through the odd/even interval method, and since it provides a full description of the 
system, the IPDF could be obtained as well. While this might work, in principle, the actual computation seems 
impractical. It would be desirable to find a more straightforward way, based on the odd/even interval method, to 
compute the IPDF. 
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APPENDIX A: STEADY STATE OF DBAP IN (N = 6)-RING 

Consider the DBAP on a ring of N = 6 sites, when there is initially an odd number of particles. The system may 
assume only one of 5 configurations, <7j (i = 1, 2, . . . , 5). These configurations and the transition rates between them 
are illustrated in Figure 4. From the figure, we see that the probabilities iTi of having state ai, obey the rate equations 

7T1 = — 7Tl + 7T2 , 
7T 2 = 7Tl — 3tT 2 + 7T 3 , 
7V3 = 27T2 — 27T3 + 27T4 , 
7T4 = 7T3 — 37T4 + 3^5 , 
7T5 = 7T4 — 37T 5 , 

where dots denote differentiation with respect to time. These equations are linearly dependent, because of the 
normalization condition X^ 71 "* ~ 1- Supplementing the rate equations with this condition, wc find the steady state: 
7Ti = 3/16, 7r 2 = 3/16, 7r 3 = 6/16, 7r 4 = 3/16, 7r 5 = 1/16. The G n may be computed for each configuration, by 
averaging over all possible locations of the n-segment (Table I). We can now compute the average G n , by weighing 
the values in Table I with the Tti found above: (G n ) — ^ i iTiG n (ai). This procedure yields the expected result, 
(G n ) = 1/2 (n = 1,2, .. . ,5). Rings of other sizes may be analyzed in much the same way, though wc found no 
obvious generalization, beyond the simple fact that (G n ) is always 1/2. 
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FIG. 1. Computation of F n (t): Empty/solid rectangles symbolize n-segments with an even/odd number of particles. 
Empty/solid circles represent empty /occupied sites. 
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FIG. 2. Change of G n with input of adjacent particle pairs: Hatched rectangles represent n-segments with initially an even 
(odd) number of particles, (a) Input of a pair inside the interval does not affect the parity, (b) Input of the pair at the edge of 
the segment, when just one particle lands inside, shifts the segment parity to odd (even). 
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FIG. 3. The crossover function f(z) of Eg. ([33[). Notice the limiting behavior for small and large z and the crossover about 
z ~ 1 evident from the plot. 
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FIG. 4. Configurations of DBAP in a (TV = 6)-ring: Empty/solid circles symbolize empty /occupied sites. The arrows and 
numbers indicate relative transition rates between configurations. 



TABLE I. G n for the configurations of DBAP in a 
(TV = 6)-ring, shown in Figure 4 
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